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PACS 42 . 65 . Tg - Optical solitons; nonlinear guided waves 

Abstract -We introduce a system based on dual-core nonlinear waveguides with the balanced gain 
and loss acting separately in the cores. The system features a "supersymmetry" when the gain and 
loss are equal to the inter-core coupling. This system admits a variety of exact solutions (we focus 
on solitons), which are subject to a specific subexponential instability. We demonstrate that the 
application of a "management" , in the form of periodic simultaneous switch of the sign of the gain, 
loss, and inter-coupling, effectively stabilizes solitons, without destroying the supersymmetry. The 
management turns the solitons into attractors, for which an attraction basin is identified. The 
initial amplitude asymmetry and phase mismatch between the components transforms the solitons 
into quasi- stable breathers. 



Introduction. Nonlinear dissipative systems subject to 
""^he parity-time iVT) symmetry feature a unique com- 
^.^bination of conservative and dissipative properties, the 
^most fundamental one being the real spectrum of eigen- 
values generated by non-Hermitian Hamiltonians 1 . In 
^pthe simplest case, such a Hamiltonian contains complex 
potentials V{x) which obeys the conjugation constraint, 
• V{x) = ]/*(— x). Recently, it was predicted and 
^^demonstrated in experiments [lO,ll that T^T-symmetric 
I systems can be implemented in optics, making use of sym- 
metric landscapes of the refractive index, combined with 
antisymmetrically placed linear amplifiers and absorbers, 
.^see also recent review 12 . The theoretical analysis of VT 
settings that may be implemented in optics was extended 



solitons. Further analysis has corroborated the existence 
of one-dimensional solitons in optical [21-^23^^25^^26^ and 



plasmonic 25 versions of the system. Moreover, stable 
fundamental solitons [27 and vortices [28", W were pre- 
dicted in similarly designed dual laser cavities with sat- 
urable nonlinearity. The advantage of this scheme in com- 
parison with the ordinary (passive) dual-core systems is 
that it may be used as the basic element of various lasing 
In particular, the PT symmetry can be 



23- 29 



^by the incorporation of nonlinearities 13 16] (in partic 



ular, of cubic gain and loss that are subject to the VT- 
invariance condition in its nonlinear form 17 18]). Soli- 



tons were predicted in some of these nonlinear systems 



19 . 



A setting which may give rise to a T^T-symmetric sys- 
tem is a Kerr-nonlinear dual-core waveguide, with the gain 
and loss separately applied to the two cores, which are lin- 
early coupled by the tunneling of electromagnetic waves 
across the separating barrier. This setting was proposed in 
Ref. 20 as a medium capable to support stable dissipative 



setups 

realized in this dual-core system by balancing the gain and 
loss in the active and passive linearly-coupled cores. This 
implies that both cores can be made of a doped lossy 
material, while the dopant atoms, if pumped externally, 
may add a linear gain to the system 24 . Then, if the 
external pump is applied to a single core, it may supply 
the gain overcompensating the intrinsic loss of this core, 
so as the net gain provides the PT balance with the other 
core, that remains lossy. The PT -invariant variant of the 
active system is especially interesting, as it corresponds 
to the lowest gain level which is capable to stabilize the 
system. 

The dual-core system makes it also possible to define 
the regime of supersymmetry^ with the strength of the 
linear coupling between the cores equal to the common 
values of the gain and loss coefficients. It is easy to find 
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a broad class of exact solutions for the dual-core system 
subject to the supersymmetry constraint, including soli- 
tons. A challenging issue is the stability of the solitons, as 
the supersymmetric system is not integrable. In this work, 
we find a specific subexponential (resonant) instability to 
which the supersymmetric solitons are subject. Then, we 
demonstrate that the solitons may be stabilized in the su- 
persymmetric system by means of the management tech- 
30-32|, which periodically switches the gain and 



nique 

loss between the cores, simultaneously reversing the sign 
of the linear coupling. The former element of the manage- 
ment scheme may be provided by the periodic switch of 
the core to which the external pump is applied, while the 
latter one may actually be implemented by periodically 
shifting the phase of one of the field components by tt, 
i.e., periodically passing the wave propagating in one core 
through TT-shifting slabs. The management scheme does 
not break the supersymmetry, making the solitons stable 
in a broad parameter region, which is identified by means 
of systematic simulations. The stabilization provided by 
this technique complies with the general fact that the man- 
agement helps stabilize solitons in a variety of waveguiding 
systems 31,32 . 

The model. The propagation of optical or plasmon 
waves in the dual-core system obeys the linearly cou- 
pled equations for the slow evolution of wave amplitudes 
u{Z,T) and v{Z^T) in the amplified and damped cores 
|20l- [251: 



iuz -\- {1/'2)utt — iju -\- Hiv = 0, 
ivz ^ {1/2)vtt ^\v\'^v ^iTv ^ K.U = 0. 



(1) 



Here Z is the propagation distance, T is the reduced time 
or transverse coordinate in the temporal or spatial system, 
respectively, 7 and F are coefficients of the linear gain and 
loss acting in the two cores, and hi is the strength of the 
linear coupling between them. The coefficients in front 
of the group-velocity dispersion (or diffraction) and Kerr 
terms in Eqs. ([T]) are scaled to be 1. The T^T-symmetry 
takes place for F = 7, while the supersymmetry condition 
is F = 7 = which casts Eqs. ([T]) into the following form: 



iuz + {1/2)utt + I'^I u — ihiu + k,v 
ivz + {1/2)vtt + \v\'^v + ihiv + Hiu 



0, 
0. 



(2) 

Obviously, substitution v {Z^T) = iu{Z^ T) transforms 
both equation ([2| into the standard NLS (nonlinear 
Schrodinger) equation, iuz + (1/2)utt + \u\^u = 0, hence 
any solution to the latter equation generates an exact so- 
lution of the supersymmetric system, the most important 
one being the fundamental soliton with amplitude 77, 

iisoi = V exp {irfZ/4) sech {r]T) , ^soi = ^^soi- (3) 

The linearization of Eqs. Q around any exact solution 
of the form v = iu leads to the following equations for 
small perturbations Su and dv: 

C {Su + iSv) = 0, £ {du - iSv) = 2m {6u + i6v) , (4) 





Fig. 1: (Color online) Instability of a fundamental soliton in the 
supersymmetric system ^ with k, = 1. The solution is gen- 
erated by the input corresponding to soliton ([3| with a small 
perturbation, uo(T) = —ivo{T) = 1.03sech(T). Panels (a) and 
(b) display the evolution of the wave components in the ampli- 
fied and damped cores, respectively (note the difference in the 
scales of the vertical axes between (a) and (b)). 



where CSu = [idz + (1/2)5tt + 2\u\'^] 6u + u'^6u'' corre- 
sponds to the linearization of the NLS equation. If u{Z^ T) 
is stable as the solution to the NLS equation (the funda- 
mental soliton deffnitely is), the ffrst equation in system 
Q produces no instability, while a resonance is antici- 
pated in the second, inhomogeneous, equation, as 5u^i5v 
is an eigenmode (corresponding to of the same operator 
C which ffgures on the right-hand side of equation. As 
follows from the commonly known linear-resonance the- 
ory, perturbations generated by the latter equation grow 
linearly in Z, rather than exponentially. Simulations of 
Eqs. Q conffrm that all the solitons are unstable, the 
character of the instability development being consistent 
with the subexponential growth, see Fig. [l] 

As said above, we aim to develop the management 
method for stabilizing solitons in the "supersymmetric" 
system, which amounts to periodically switching the gain 
and loss between the cores, with the simultaneous rever- 
sal of the sign of the coupling coefficient. The accordingly 
modified equations ([2| are 



iuz + {1/2)UTT + 1^1^^ 
+/T:sgn {sin (27rZ/L)} {—iu + 'i;) = 0, 

ivz + {1/2)vtt + \v?v 
-^nsgii {sin (27rZ/L)} {iv + ii) = 0, 



(5) 



where L is the management period. We stress that the 
management does not break the supersymmetry of the sys- 
tem, and, accordingly, wave forms with v{Z^ T) = iu{T, v)^ 
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where u is any solution to the usual NLS equation, remain 
exact solutions of system ([5|. 

The stabilizing effect of the management can be ex- 
plained by its action on perturbations around the exact 
solutions. In particular, the replacement of constant Hi in 
the second equation in system Q by the periodically flip- 
ping coefficient detunes the destabilizing resonance. The 
detuning is small for a large management period, hence 
the stabilization effect is expected to be weak in this case, 
which is indeed observed, see Fig. [3] below. Previously, 
effects of the periodic modulation of the coupling strength 
on the stability of solitons were studied in the model of 
the usual (conservative) dual-core systems 30,32 . Those 
effects may be rather involved, but, generally speaking, 
they do not stabilize the symmetric solitons against the 
symmetry-breaking bifurcation, as rapid modulations ef- 
fectively weaken the linear coupling between the cores, 
thus making the symmetric solitons less robust against the 
spontaneous symmetry breaking. Thus, the effects of the 
management mechanism acting on the "supersymmetric" 
solitons are different from the earlier studied ones. 

The stabilization of solitons by means of the manage- 
ment. Simulations of Eqs. ([5| reveal that the periodic 
application of the switching not only stabilizes the exact 
soliton solutions against small perturbations, but actually 
turns the solitons into strong attractors for inputs signif- 
icantly different from the exact solutions. Typical exam- 
ples illustrating the rapid convergence of the inputs, taken 
as 

uo{T) = A sech{T),vo = iuo, (6) 

with amplitudes A differing by +20% or —10% from the 
value = 1 in the respective exact soliton solution ([3|, 
into solitons are displayed in Fig. [2] Note that the prop- 
agation distance shown here in units of Z actually coin- 
cides, within the order of magnitude, with that measured 
in units of the soliton's dispersion length, therefore the 
propagation distance displayed in the figure is extremely 
large, making the full stabilization of the solitons obvious. 

While the stabilization seems rather obvious in the case 
of the short-period management (as in that case the fre- 
quent switch between the loss and gain may nearly average 
them to zero), the model provides conspicuous stabiliza- 
tion even for very large management periods. The results 
of systematic simulations are summarized in Fig. [3j which 
displays stability domains (areas between the curves) in 
the plane of the management period L and input ampli- 
tude A (defined as per Eq. (|6|), for different values of the 
coupling parameter k in Eqs. ([5|. While it is natural that 
the stability area expands with the decrease of (which 
corresponds to weaker perturbations), it is worthy to note 
that the stability regions remain appreciable even for very 
strong coupling {hz = 10). Note that (as was mentioned 
above) the stability areas may be realized as attraction 
basins for the stable solitons which self-trap from input 
pulses (|6|. 

Amplitude B of the self-trapped soliton, that has the 
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Fig. 2: (Color online) The self-trapping of inputs (|6| with A = 
1.2 (a) and A = 0.9 (b) into stable solitons, demonstrated by 
simulations of Eqs. ^ for k = 1 and L — 2. The evolution of 
the t'-component is similar to that displayed here for u{Z,T). 




Fig. 3: (Color online) The soliton stability regions in the (L, A) 
plane for different values of coupling parameter k, (see Eqs. ^ 
and ([6|). The inset shows amplitude B of the emerging stable 
soliton as a function of input amplitude A. 
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Fig. 4: (Color online) The evolution of the peak powers of 
both field components in the case of the initial conditions with 
asymmetric amplitudes, uq = 1.05sech(T), vq = 0.95i sech(T). 
The coupling parameter is At = 0.25. (a) The system under the 
action of the management with period L — 2. (b) The system 
without the management. 

form of \u{T)\ = \v{T)\ = B sech(5T), is a function of 
the input's amplitude A of initial condition ([6|, as shown 
in the inset of Fig. |3] It is worthy to note that, unlike the 
stability area, dependence B{A) is weakly affected by the 
variation of L and hi. 

Obvious peculiarities of the exact supersymmetric so- 
lutions, with V = iu^ are the equality of the amplitudes 
of the two components, and the fixed phase shift between 
them, 7r/2. The use of inputs which break either of these 
relations impedes the self-trapping into exact solitons. In 
such a case, the management gives rise to overall-stable 
two-component solitary breathers^ which feature persistent 
irregular oscillations and extremely slow evolution of the 
average peak power. A typical evolution history produced 
by the input with the asymmetry in amplitudes is illus- 
trated by Fig. Qa). For comparison. Fig. 4(b) shows 
that the same input pulse quickly collapses in the system 
without management. 

The solitons in the supersymmetric system are also sen- 
sitive to a mismatch in the phase shift between the two 
components. The simulations with input 

uo{T) = sech(T), ^0 = iexp(z(/))^io(T), (7) 
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Fig. 5: (Color online) The same as in Fig. [4] but for k, — 
1 and initial condition ([7| with phase mismatch between the 
components: (j) — 0.257r in (a) and (j) — O.Stt in (b). 

where (j) is the phase mismatch (cf. Eq. (|3|), generate 
quasi-stable two-component breathers, which remain cou- 
pled over a finite but long propagation distance, if (j) is not 
too large. An example for (j) = 0.257r is displayed in Fig. 
[sja). On the other hand, the increase of (j) in the same case 
to (j) = O.Stt leads to quick decoupling of the components 
and destruction of the two-component soliton, as shown 
in Fig. [lb). 

Collecting results of the simulations, it is possible to 
identify a critical value of the mismatch, ^cr, which sepa- 
rates the quasi-stable transmission of the two-component 
soliton and its quick splitting. Fixing = 1, is dis- 
played as a function of the management period, L, in Fig. 
[6] It is observed that the phase mismatch produces the 
strong destabilization effect at small and large values of L, 
while a flat optimum region is identified at 1 ^ L ^ 2.5. 
For 1/ < 1, detailed consideration of the numerical results 
explains the destabilization as follows: small L affects, first 
of all, the coupling between the components, practically 
nullifying it due to the averaging over frequent flips of the 
sign of hi. Thus, the system splits into almost decoupled 
components, and the one which starts its evolution under 
the action of the loss cannot recover after the initial de- 
cay, despite the subsequent switch of the loss into the gain 
(similar to what is seen in Fig. jsjb)). In the opposite 
case of large L, the splitting of the components is sim- 
ply explained by the fact that one of them suffers strong 
attenuation during the first period of the management. 

Conclusion. We have introduced a system based on 
the dual-core coupler with the intrinsic Kerr nonlinearity, 
which features the super symmetry, as an extended ver- 
sion of the VT symmetry, with the gain and loss in the 
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Fig. 6: (Color online) The critical value of the initial phase mis- 
match, (/)cr, for the input ([7|, versus the management period, 
L, for At = 1 in Eq. ([5f. 

two cores equal to the strength of the inter-core coupling. 
The supersymmetry allows one to find a broad class of ex- 
act solutions, including solitons, for which the stability is 
a crucial issue. The exact solutions are vulnerable to the 
specific subexponential instability, which destroys all the 
solitons. We have demonstrated that the management, 
defined as periodic synchronous switch of the signs of the 
gain, loss, and linear coupling readily stabilizes the soli- 
tons, keeping the system's supersymmetry intact. Actu- 
ally, the management makes solitons attractors, for which 
we have identified the attraction basins. The amplitude 
asymmetry and phase mismatch between the two compo- 
nents in the input pulse turns the stationary solitons into 
quasi-stable breathers. We have identified the stability 
area for the breather with respect to the phase mismatch 
and management period. 

As an extension of this study, it may be relevant to 
consider interactions between the stabilized solitons. A 
challenging problem is to analyze two-dimensional double- 
core supersymmetric systems. 
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